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Abstract 

The MIEZE (Modulation of Intensity with Zero Effort) technique is a variant of neutron resonance spin echo (NRSE), which has 

proven to be a unique neutron scattering technique for measuring with high energy resolution in magnetic fields. Its limitations 

in terms of flight path differences have already been investigated analytically for neutron beams with vanishing divergence. In the 

C^ present work Monte-Carlo simulations for quasi-elastic MIEZE experiments taking into account beam divergence as well as the 

^Sl sample dimensions are presented. One application of the MIEZE technique could be a dedicated NRSE-MIEZE instrument at the 

5_i European Spallation Source (ESS) in Sweden. The optimisation of a particular design based on Montel mirror optics with the help 

Qh of Monte Carlo simulations will be discussed here in detail. 
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^_^ 1. Introduction 

-(— > 

^ One major field of current solid state physics is complex 
I magnetic systems. Examples are the structure and dynamics 
^ of topological spin structures like skyrmions in metals, semi- 
.^^ conductors and ferroelectric insulators [1, 2, 3, 4]. Others are 
y^ the dynamics of spin waves in ferromagnets like EuO and Fe 
O [5]. Further interest focuses on the slow magnetic dynamics of 
cZ) magnetic monopoles and Dirac strings in spin ice [6]. Gener- 
K^ally, excitations with lifetimes from the picosecond to the mi- 
'^ crosecond range near quantum phase transitions are of utmost 
L_J interest for the understanding of the magnetic dynamics in com- 
plex materials, especially with regard to possible applications in 
^~^ spintronics. 

_. The standard technique to measure quasi-elastic dynamics 
.^ in the higher ps and ns time range is Neutron Spin Echo (NSE) 
(^ [7]. As a Larmor precession technique it is sensitive to mag- 
l/~J netic fields and can only be used with magnetic fields at the 
T^ sample position by applying rather complicated spin manip- 
f^ ulation resulting in a loss of intensity [5]. The complemen- 
C^ tary Neutron Resonance Spin Echo (NRSE) method invented 
". . by Golub and Gahler [8, 9] shares the sensitivity to magnetic 
^ fields or depolarising samples, but can be evolved straightfor- 
k> wardly into the MIEZE technique [9, 10, 11], which is then 
%_J independent of sample depolarisation effects, as all spin manip- 
d ulations are performed before the sample. 

As MIEZE is effectively a time-of-flight method it is sen- 
sitive to differences in the length of the neutron flight path. 
This drawback has been studied analytically in reference [12] 
for beams with zero divergence. But for analysing MIEZE data 
and designing MIEZE instruments the effect of divergence and 
the influence of the neutron path on the resolution is required. 



In this study we will show that using the Monte Carlo simu- 
lation package McStas [13, 14] together with specifically writ- 
ten components for the MIEZE setup, both earlier analytical 
and experimental results can be reproduced and furthermore be 
used to assess the suitability of more elaborate instrument de- 
signs. MIEZE is an excellent technique for high energy res- 
olution measurements: For small samples spin-echo times of 
the order of a microsecond are possible with currently available 
components such as RF spin flipper coils and TOF detectors. 

2. MIEZE technique 

As a variant of NRSE, the MIEZE technique omits both tt- 
flipper coils after the sample [9]. The coils before the sample 
are separated by a distance Li and are driven at two different 
frequencies, wi < W2. After passing the RF-flippers, which are 
working in resonant 7r-flip mode without bootstrap, the neutron 
spin phase depends on the flight time and the neutron velocity 
V. At a distance L2 given by 



wi ■ Li = (co2 - wi) ■ L2 



(1) 
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the spin modulation of all velocity groups adds up in phase. 
Eq. (1) is called the MIEZE condition. The key point of MIEZE 
is now that the spin modulation can be converted to an inten- 
sity modulation by placing a polariser anywhere after the sec- 
ond coil, in particular before the sample. Then a time-resolved 
detector is required, which registers a sinusoidally modulated 
signal of the frequency um I{t) oc cos,{(±iMt), where com - 
2(w2 - cjDi) is the MIEZE frequency. 

If the phase at the detector is not uniform, which can be 
caused by inelastic scattering at the sample or by path length 
differences after the polariser, the modulation is damped and 
the signal can be described by /(f) - B + A cos(wmO- For neg- 
ligible path length differences the MIEZE contrast C - AjB 
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is a measure of the inelasticity of the scattering, which is the 
MIEZE equivalent of the measured polarisation in NSE/NRSE. 
It can be shown [15] that with an ideal setup, the contrast is 
the cosine Fourier transform of the scattering law S (q, co) of the 
sample 
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C{q,TM) - I d(xiS{q,to) cos(wtm), 



(2) 



with the sample placed at the position Ls before the detector. 
With m being the mass of the neutron, the MIEZE time tm is 
given by 



Tm 



TudmLs 



mv^ 



(3) 



This is the Fourier time of the instrument, equivalent to the spin- 
echo time of NSE/NRSE [15]. 

For a real MIEZE instrument, the contrast of the modula- 
tion is reduced by imperfections of polarisers, coils, detector 
and by path length differences. The reduction in contrast C due 
to different sample shapes and spin-echo times was studied the- 
oretically in reference [12]. To benchmark the quality of the 
Monte-Carlo simulation of MIEZE, we reproduce and verify 
the previous work in the next section. 

3. MIEZE and Monte-Carlo simulations 

3.1. Reproduction of previous results 

For the reproduction of the previous results from the refer- 
ence [12] we used the Monte-Carlo simulation package McStas 
[13, 14] simulating a standard MIEZE setup without any addi- 
tional neutron optics and the sample only. 

The contrast of the MIEZE signal is reduced by several ef- 
fects of which we discuss the most important ones, namely the 
thickness of the absorber in the detectors, the resolution of the 
detector and the geometry of the sample in more detail. The 
effect of imperfections in the NRSE coils on the spin and the 
dependence on the MIEZE frequency um is not taken into ac- 
count here as the coils are simulated as ideal RF spin flippers. 

The thickness Ad of the detector modifies the registered 
contrast by introducing a phase shift in equation (2) 
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Mj Jo 



dx S {(jj) cos {(jjtm + cl>mx/v) . (4) 



A phase shift of In over the whole thickness of the detector 
would completely destroy the contrast, whereas a phase shift of 
n/2 would not reduce it by more than 10 %. For example in 
the CASCADE detector [16, 17, 18] we used in our MIEZE ex- 
periments, where the thickness of the detection planes is Ad = 
2/im, a 10 % loss of contrast happens at MIEZE frequencies 
of wm = 621 MHz, COM = 311 MHz and com = 207 MHz for 
A - 5 A, A - lOA and A - 15 A, respectively. For the re- 
spective wavelengths this corresponds to a maximum possible 
MIEZE time of t^ = 1 yus, tm = 4/is and tm = 9//s, which is 
far above the values considered here (see Figures 6 and 7) thus 
justifying the assumption of a thin detector in the simulations. 

To counter the effect of contrast reduction due to path length 
differences within the plane of the detector especially at higher 



MIEZE times tm an area of 1 mm^ is used in the simulation. 
This is in the order of the spatial resolution of a CASCADE 
detector. 

As the thickness and spatial resolution of the detector there- 
fore do not contribute significantly to the contrast reduction, we 
only have to consider the reduction factor due to the sample in 
the following. 

In the simulations the contrast is normalised to its value at 
(7 = 0. Vanadium is used as a purely elastic scatterer. Fig. 1 
shows a comparison of the simulated and calculated [12] re- 
duction factors due to different sample shapes in fixed and re- 
flecting geometry (where reflecting is defined as a rotation - 
0/2 - 90° in the direction of 6, see Fig. 1) For the spherical 
and cylindrical samples absorption is neglected. This simplifi- 
cation is necessary because the analytical formulae in [12] do 
not take into account the additional contrast reduction caused 
by inhomogeneous absorption due to sample shape. The results 
demonstrate an almost perfect agreement between the analyti- 
cal and the simulated results. 

Fig. 2 shows the simulated (solid lines), calculated (dashed 
lines) and measured (symbols with error bar) reduction factors 
for a cuboid sample in fixed geometry for different spin-echo 
times Tm- There is an excellent agreement between simulation 
(this work), calculation and measurement (from [12] ). 

Summarising, the effects of contrast reduction due to the ge- 
ometry of the sample calculated using an analytical method are 
very well reproduced by the McStas simulations using a stan- 
dard MIEZE setup with no additional neutron optics (compare 
Figs. 1 and 2 in this paper to Figs. 6 and 7 in [12]). 

3.2. Compensation for detector misalignment 

In a real experiment with a total length in the order of sev- 
eral tens of metres a detector misalignment relative to the coils 
violating the MIEZE condition, will lead to a reduction of MIEZE 
contrast. We briefly show that this effect can easily be compen- 
sated by adjusting one of the RF frequencies. 

A detector that is misaligned by an offset of Al away from 
L2 does not fulfill the MIEZE condition (Eq. 1) anymore. How- 
ever, the signal can be recovered by adjusting the frequency of 
one of the 7r-flipper coils. This can immediately be seen if we 
substitute L2 by L2 + ^l and wi by coi + Awi in eq. (1) and solve 
for Aw 1: 



Awi — (co2 — coi) 



L\ ■Cl>2 

Al + Li + L2 



(5) 



According to Eq. (5), the first spin flipper frequency io\ has 
to be adjusted by Aco\ to compensate for a detector misalign- 
ment of Al. The simulations (Fig. 3) show that the signal can be 
perfectly recovered. 

4. A dedicated NRSE-MIEZE instrument at the ESS 

In the following we show how the newly developed simula- 
tion tools can be used to demonstrate the feasibility of a focus- 
ing MIEZE beam line for various q and tm parameters. 
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Figure 1 : Top: The different sample shapes are shown. In the reflecting geome- 
try the sample is rotated at an angle = 8/2 - 90° with respect to the scattering 
angle as indicated in the figure. Bottom: Comparison of the sample reduction 
factor for different geometries with constant sample volume V = 100 mm' us- 
ing the analytical results from [12] (dashed lines) and simulated results (solid 
lines). Parameters are as given in [12]: oim = 2n ■ 200kHz, corresponding to 
tm = 1.15 ns, Li = 1 m, L2 =2m, L, = 0.8m, .i = 10.4 A. 



We will propose a combined MIEZE/NRSE instrument for 
measuring slow dynamics, in particular from depolarising sam- 
ples and samples in extreme environments (magnetic fields, pres- 
sure) at the planned European Spallation Source (ESS) in Lund, 
Sweden. This source will provide neutron pulses with a repe- 
tition rate of 14 Hz, a pulse length of 2.8 ms and a total power 
of 5 MW. The peak intensity will reach 30 times the average 
intensity of the ILL. A sketch of the planned instrument config- 
uration is shown in Fig. 4. It consists of a beam extraction and 
guide system based on Montel optics [19, 20, 21, 22], which al- 
lows a decoupling of the selection of the size and the divergence 
of the incident beam. The distance between the focal points of 
each Montel mirror is projected to be Li = 16 m. The reflecting 
mirrors are assumed to be 9.6 m long. The goal is to trans- 
port only those neutrons to the sample that will be useful for 
scattering from the sample, i.e. the beam size is equivalent to 
the sample size and the divergence is such that it is compatible 
with the requirements of the experiment. All other neutrons are 
absorbed before, therefore the background will be minimised 
at the sample position and the low radiation load on the trans- 
port system for the neutrons will yield a longer life time of the 
components. 

The chopper system will be designed such that it provides 
neutrons for wavelength bands 3 < /I < 9.5 A and 9 < A < 15.5 
A by changing the phase relation between the two choppers. 
They are located right after the biological shielding of the ESS 
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Figure 2: Comparison of the sample reduction factor for a cuboid in fixed ge- 
ometry showing the reduction factor for the simulation results (solid lines), the 
analytical results from [12] (dashed lines) and experimental (symbols with eiTor 
bars) results from the instrument MIRA also performed by [12]. Parameters are 
as given in [12]: sample width 25 mm, thickness 5 mm, A = 10.4 A, Ci>m ranges 
from 2n ■ 46 kHz to In ■ 200 kHz, yielding a MIEZE time tm from 0.26 ns to 
1.15 ns, respectively. Instrument lengths are: Li = 1 m, L2 = 2 m, L, = 0.8 m. 



(6 m downstream the moderator) and in front of the sample. 
Therefore, neutrons spanning the whole wavelength band with- 
out frame overlap and using all pulses from the moderator are 
used. The lower wavelength band around 6.25 A defines the to- 
tal length of the instrument (Liot - 45 m) from the moderator to 
the sample position, which is located 26 m after the biological 
shielding (Fig. 4). 

The coil system consists of two NRSE coils: one in front of 
the sample and one located between the two Montel mirrors of 
the guide system. The 16 m long precession region is shielded 
by y[/-metal for providing a zero field region. The first Mon- 
tel mirror polarises the beam. In front of the sample will be 
a removable spin analyser. Running the two coils at different 
RF-frequencies one obtains a time modulated intensity at the 
detector Two options are possible: 

i. For the small angle scattering regime, a detector of 1 m^ 
with a spatial resolution of 2 mm will be placed 10 m be- 
hind the sample allowing a time resolution of the order of 
yus (see Figs. 6 and 7). Such a detector can be realised for 
example by combining position-sensitive detectors such as 
the CASCADE system [16, 17, 18]. It also conserves the 
MIEZE contrast as described above, 
ii. For larger values of q, a set of NRSE coils behind the sam- 
ple will correct for the path length differences in the beam. 
Using an analyser bank over an angular range of 90° and a 
detector distance of 2 m (with a relaxed resolution of 1 cm), 
this setup will allow for a (7-range up to 3 A"' with a time 
resolution exceeding 100 ns. Note that this configuration 
has not yet been simulated as it needs a different simula- 
tion method. 

The Monte-Carlo simulations as described in section 3.1 
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Figure 3: Compensation of a detector offset of 2 cm by adjustment of the first 
spin flipper frequency. The left curve shows the signal for an aligned detector, 
the right curve demonstrates that the signal can be recovered for a misaligned 
detector by changing the frequency /m = (jjMl2n of the second coil according 
to Eq. (5). Instrument parameters are: L| = 16m, Li = 16m, L, = 13m, 
A=\5k. 




Figure 4: Schematic view of the MIEZE instrument using Montel optics. Mi 
and Ml are reflecting mirrors with a length of 9.6 m each, which is 60% of 
the distance L\ between the n'-flipper coils Ca and Cg, A and S indicate the 
analyser and the sample, respectively. The instrument parameters used in the 
simulations are: Distance from the virtual source with an area of 2 x 2 mm^ 
to the coil Ca'- Lq = 13 m, distance L[ = 16 m, distance from coil Cb to the 
detector: L2 = 16 m, distance from coil Cb to the sample: Lc = 3 m, distance 
from the sample to the detector: L , = 13 m. 



were now used to simulate the signal for the small angle regime 
of the proposed instrument i.e. option (i) above. The simulation 
consists of the ESS source, the instrument including the Montel 
optics, the MIEZE coils and the detector Perfect MIEZE coils, 
Montel mirrors and no gravity were assumed. Note, that the to- 
tal contrast reduction is obtained by multiplying the individual 
contrast reductions due to the focusing optics (Fig. 6)) and the 
sample shape (Fig. 7). 

An optimisation of the inclination of the NRSE coils due 
to the chosen mirror configuration (see Fig. 5) was performed. 
The influence of the flight path differences on the resolution 
was significantly reduced for an optimal inclination of the coils 
by 1.65 degrees. A misplacement of the detectors would degen- 
erate the signal, but can be compensated by the adjustment of 
the MIEZE frequency as demonstrated in Fig. 3. Using a wave- 
length of 15 A and the optimal MIEZE frequency, which allow 
for maximal MIEZE times of up to a microsecond, the simula- 
tion results show no relevant decrease of the MIEZE contrast at 



Figure 5: Dependence of the contrast on the tilt of the n'-flipper coils. The 
optimal coil rotation of 1.65" is independent of A. 



the optimal inclination of the coils (see Fig. 6). Taking into 
account the reduction factor due to the sample shape the acces- 
sible q-TM space for two typical sample geometries used in the 
field of magnetism is plotted in Fig. 7. The results demonstrate 
that the large time resolution up to one microsecond of the in- 
strument can be preserved for selected sample geometries. 

In conclusion, we succeeded to prove that Monte-Carlo sim- 
ulations are a powerful means for the optimisation of spectrom- 
eters using the MIEZE technique. Choosing a favourable sam- 
ple geometry can preserve the excellent time resolution of the 
instrument. Furthermore, the use of focusing techniques and 
Montel optics will lead to a low background and a tremendously 
reduced radiation load on the beam components. Moreover, no 
mirror optics and no choppers are required inside the biological 
shielding thus facilitating the maintenance of a MIEZE beam 
line. Adjusting the MIEZE frequency will compensate for pos- 
sible detector misplacements. We have shown that path length 
differences can be effectively compensated by properly tilting 
the ;7r-flippers in a focussing NRSE/MIEZE instrument. 

We have demonstrated that a powerful, low background and 
low cost MIEZE/NRSE beam line at the ESS can be realised us- 
ing beam components that are already available now and are in 
regular use at the beam line RESEDA at FRM II [23] . The time 
resolution of a MIEZE/NRSE at large angles will be of the same 
order as for NSE SPAN-like instruments [24]. However, going 
to small q will add the benefit of performing experiments from 
depolarising samples and in magnetic fields with a significantly 
improved time resolution. 
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Figure 6: Contrast vs. MIEZE time tm for various wavelengtlis for the small 
angle regime without inclusion of the effects of sample shape and gravitation. 
The instrument parameters are; L| = 16m,L2 = 16m, Lj = 13m. The coating 
of the Montel optics is m = 5 times the critical angle of natural nickel. 



Figure 7: Time resolution limitations from the sample shape, not including 
effects of the focussing mirrors. The R = 50% contour for a cuboid sample in 
reflecting geometry with a thickness of 2 mm (typically used for soft mater) and 
a cylindrical sample with a radius of 3 mm (typically used for hard condensed 
matter and for experiments under extreme conditions) are shown for difterent 
wavelengths. Instrument parameters are: Li = 16m, L2 = I6m, L, = 13m. 
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